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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 

RESEARCH  MEMORANDUM 

APPLICATION  OP  STATISTICAL  THEORY  TO  BEAM-RIDER 
GUIDANCE  IN  THE  PRESENCE  OF  NOISE. 

I - WIENER  FILTER  THEORY 
By  Elvood  C.  Stewart 

SUMMARY 


A study  has  been  made  of  the  application  of  Wiener  filter  theory  to 
the  design  of  a beam-rider  guidance  system  operating  in  the  presence  of 
glint  noise.  Target  and  missile  motions  are  restricted  to  the  same  plane. 
The  Wiener  theory  is  then  used  to  establish  the  theoretical  lover  limit 
of  root-mean- square  error  and  the  corresponding  desired  transfer-function 
characteristics.  It  is  shown  that  although  the  practical  achievement  of 
these  results  is  restricted  by  saturation  effects,  the  theory  is  useful, 
with  suitable  modifications,  as  a guide  in  system  design.  Such  modifica- 
tions have  been  applied  to  the  design  of  Bystems  for  which  the  optimum 
filtering  is  placed  either  in  the  missile-control  system  or  in  the  tracking 
radar.  The  error  performance  of  these  systems  for  different  noise  magni- 
tudes is  presented.  Other  considerations  such  as  servo  energy  requirements 
are  briefly  discussed. 


INTRODUCTION 


In  the  design  of  a missile-guidance  system  certain  standard  criteria, 
such  as  the  fastest  response,  are  not  always  the  most  useful.  This  is 
particularly  true  in  the  case  of  guidance  systems  which  are  forced  to 
operate  in  the  presence  of  certain  random  unwanted  disturbances  known  as 
noise.  The  reason  for  this  is  that  the  effects  of  the  noise  can  seriously 
reduce  the  probability  that  the  missile  will  hit  the  target.  Furthermore, 
the  sources  of  noise,  being  dependent  on  the  target  characteristics,  can- 
not be  eliminated.  Consequently  the  guidance  system  should  be  designed 
to  minimize  the  miss  distance  even  when  the  noise  is  present.  This  prob- 
lem will  be  considered  here. 

There  are  two  possible  design  approaches.  In  the  first,  the  form 
of  the  system'  (i.e.,  the  transfer  function)  is  assumed  and  an  attempt  is 
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made  to  adjust  the  existing  parameters  so  as  to  reduce  the  effects  of 
noise  (ref.  l) . This  procedure  is  not  only  difficult  to  apply,  but  the 
ultimate  performance  is  limited  by  the  assumed  form.  In  the  second 
approach,  the  form  as  well  as  the  parameters  of  the  transfer  function 
are  determined  so  that  the  noise  has  the  least  possible  effect  on  the 
performance  of  the  missile;  an  attempt  is  then  made  to  devise  a system 
which  has  a transfer  function  approximating  this  optimum  transfer  func- 
tion. The  latter  approach  will  be  adopted  here. 

The  problem  of  determining  optimum  transfer  functions  has  been 
encountered  previously  in  the  communication  field,  and  as  a result  of 
this  encounter  a statistical  theory  known  as  Wiener  filter  theory  has 
been  developed.  By  the  use  of  this  theory  it  is  possible  to  determine  a 
unique  optimum  transfer  function  which  will  result  in  a theoretical  lower 
limit  of  mean-square  error  between  the  desired  and  the  actual  missile 
position.  Very  little  work  has  been  done  in  the  application  of  this 
theory  to  beam-rider  guidance.  Previous  works,  references  2 and  3,  have 
been  confined  to  simple  homing  systems.  The  purpose  of  the  present  paper, 
however,  will  be  to  investigate  the  applicability  of  this  theory  to  a 
beam- rider  guidance  system. 

In  the  application  of  this  theory  to  missile  guidance  it  is  necessary 
to  make  certain  assumptions.  Foremost  of  these  is  the  assumption  that  the 
target  and  missile  move  in  the  same  plane , taken  in  this  report  to  be 
horizontal.  Other  assumptions,  such  as  those  relating  to  the  class  of 
target  maneuvers  and  noise,  are  discussed  in  the  text.  Within  these 
restrictions,  however,  the  theory  may  be  used  to  obtain  a measure  both 
of  the  error  performance  that  might  be  expected  and  the  difficulties  to 
be  overcome  in  order  to  realize  this  performance. 


SYMBOLS 

N noise  magnitude  or  zero  frequency  spectral  density,  ft2/radian/sec 
TN  time  constant  of  the  noise  spectrum  shaping  filter,  sec 
Y0  optimum  closed- loop  transfer  function 

a acceleration  of  target  maneuver,  ft/sec2 

k twice  the  average  switching  rate  of  target  acceleration,  1/sec 

y apparent  target  displacement  from  true  target  center  due  to  noise, 
H ft 

y„  missile  displacement  from  a space  reference,  ft 

M 
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target  displacement  from  a space  reference,  ft 

spectral  density  of  noise  displacement  ft2/radian/sec 

spectral  density  of  target  displacement  y^,  ft2/radlan/sec 

error  between  target  and  missile  position,  y,p  - y^,  ft 

component  of  error  e due  to  noise,  ft 

component  of  error  e due  to  target  motion,  ft 

control- surface  deflection,  radians 

angle  of  yaw,  ^ i ms 

optimum  open-  transfer  function 

open- loop  transfer  function  of  system  approximation  to  n0 
angular  frequency,  radians/ sec 


GENERAL  CONSIDERATIONS 


* Of  the  many  sources  of  noise  which  may  exist  in  a guidance  system 

utilizing  radar  detection,  glint  noise  is  one  of  the  most  serious.  Glint 
noise  is  a term  that  is  used  to  describe  a shift  in  the  apparent  center 
- of  a target  as  determined  by  a tracking  radar.  It  Is  due  basically  to 

the  variable  reflection  characteristics  of  aircraft  targets  and  arises 
from  the  relative  movement  of  the  various  reflecting  surfaces.  Since  the 
radar  utilizes  the  reflected  signal  to  determine  target  position,  vari- 
ations in  the  reflected  signal  are  interpreted  by  the  radar  as  a shift 
f in  the  target  center.  This  type  of  noise  is  particularly  Important  since 

It  is  due  fundamentally  to  the  target  characteristics  and  therefore  cannot 
i be  eliminated  by  any  known  radar  improvements.  The  situation  is  illus- 

trated in  figure  1(a)  where  the  true  target  position  is  indicated  as  y_, 
and  the  glint  noise  is  represented  by  the  displacement  yN. 

I The  present  report  Is  restricted,  for  the  sake  of  simplicity,  to  a 

; two-dimensional  study  in  which  the  target  and  missile  move  in  a horizontal 

i plane.1  The  guidance  system  is  considered  to  be  of  the  beam-rider  type, 

illustrated  in  figures  1(a)  and  l(b).  It  should  be  noted  here  that  dis- 
placements are  referred  to  a fixed  space  reference.  The  function  of  the 
guidance  system  is  to  make  the  missile  position  y^  coincide  as  closely 
1The  complete  three-dimensional  problem  would  require  a more  complex 
\ . analysis  than  used  herein.  Possibly  either  the  present  theory  or  Wiener's 

| theory  for  multiple  time  series  (ref.  4)  could  be  applied  to  this  case. 
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as  possible  with  the  actual  target  position  y^,.  This  requirement  would 
he  relatively  easy  to  satisfy  if  the  tracking  radar  could  locate  the 
target  precisely.  However,  because  of  the  radar  noise  the  only  informa- 
tion available  to  the  guidance  system  is  the  apparent  target  position  as 
illustrated  in  figure  1(b) . For  this  reason  the  task  is  much  more  diffi- 
cult and  the  missile  position  may  deviate  considerably  from  the  actual 
target  position.  The  difference,  y™  - y»,  is  denoted  on  this  figure  by 
the  error  €,  which  should  be  minimized  in  sane  manner. 

A suitable  criterion  for  Judging  system  performance  depends  primarily 
on  the  manner  in  which  the  system  operates.  In  the  case  of  the  beam-rider 
system  the  missile-to-target  range  is  not  normally  transmitted  to  the 
missile  so  that  the  missile  cannot  know  when  the  target  will  be  reached. 
Hence  the  error  should  be  minimized  for  all  values  of  range  or,  equiva- 
lently, of  time.  A mathematically  convenient  criterion  which  does  not 
involve  weighting  with  respect  to  missile  travel  time  is  the  mean-square 
time  average  of  the  error.  This  criterion  will  be  used  herein. 

Th-;  design  of  a system  normally  depends  on  the  inputs  to  be  encoun- 
tered, in  this  case  the  target  motion  and  noise.  Because  of  their  random 
nature  it  is  not  convenient  to  define  these  quantities  explicitly  as 
functions  of  time,  and  statistical  descriptions  are  more  suitable.  Since 
it  is  generally  believed  that  the  target  motion  and  noise  are  uncorrelated 
they  will  be  described  independently.  What  follows  is  a brief  discussion 
of  these  inputs. 

Intensive  effort  has  been  devoted  in  recent  years  to  the  measurement 
of  radar  glint  noise.  References  3 and  5 through  8 are  typical  of  such 
work.  The  quantity  of  most  general  interest  in  these  measurements  is  the 
displacement  of  the  apparent  center  of  the  target  from  the  time  center, 
or  yN,  shown  in  figure  1(a) . This  quantity  can  be  defined  statistically 
by  means  of  (l)  the  amplitude  distribution  and  (2)  the  power  spectral  den- 
sity. Although  the  determination  of  these  quantities  is  somewhat  uncertain, 
it  is  generally  found  that  the  amplitude  distribution  is  approximately 
Gaussian  and  that  the  spectral  density  can  be  adequately  represented  by 


N 

Tj|2U)2  + 1 


(1) 


Spectra  obtained  from  any  one  individual  experiment  may  deviate  somewhat 
from  this  form  but  it  is  generally  considered  that  the  above  character- 
istic represents  a reasonable  average  of  many  different  experiments. 
Examination  of  noise  spectra  indicates  that  the  break  point  (l/2*TN)  is 
on  the  order  of  several  cycles  per  second  and  for  this  report  will  be 
taken  to  be  6 cps  corresponding  to  TN  =°  .0265  second.  The  magnitude 
of  the  spectrum,  N,  depends  on  factors  such  as  target  size  and  target 
aspect  so  that  the  guidance  system  is  usually  forced  to  operate  over  a 
wide  range  of  magnitudes..  This  range  may  extend  from  7 ft2 /radian/sec 
for  smal  1 targets  up  to  around  30  ft2/radian/sec  for  large  bombers. 
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Because  of  the  variation  in  magnitude  due  to  uncertainty  in  target  size 
and  aspect  as  veil  as  uncertainty  in  the  noise  measurements,  this  factor 
becomes  of  real  and  practical  importance  and  will  be  considered  herein. 

In  considering  target  maneuvers  it  is  difficult  to  say  exactly  how 
a target  will  maneuver  when  under  attack.  However,  a reasonable  situa- 
tion might  be  one  in  which  the  target  is  merely  aware  of  the  attack  and 
therefore  maneuvers  in  some  severe  manner  to  avoid  being  hit.  Here  it 
will  be  assumed  that  the  target  maneuvers  laterally  with  maximum  acceler- 
ation alternately  in  opposite  directions.  The  duration  of  each  acceler- 
ation will  be  a random  function  determined  by  some  distribution.  A 
reasonable  distribution  which  leads  to  an  easily  handled  spectral  density 
is  the  Poisson  distribution  [ (l/3?)exp(-T/$) ] where  T represents  the  time 
and  T the  average  time  between  switching  of  the  acceleration.  As  shown 
in  reference  9 the  spectral  density  of  the  target  acceleration  is  then 
described  by 


$ _ ka2 

&T  + k2) 

Here  the  quantity  a represents  the  magnitude  of  the  target  acceleration 
normal  to  the  beam,  and  k is  twice  the  average  switching  rate,  k = 2/$. 
The  spectral  density  of  the  target  displacement  is  then  given  by 

•T  a ®yT  = „W4(^  + k2)  (2) 

It  would  appear  that  there  is  a problem  here  concerning  the  existence  of 
this  spectral  density  because  of  the  u4  in  the  denominator.  However, 
it  can  be  shown  that  it  is  possible  to  use  this  representation  for  pur- 
poses of  computations  (see  Appendix  A) . For  the  tail-chase  maneuver  to 
be  used  in  a later  example  the  target  acceleration  is  specified  to  be 
±1  g at  an  average  period  of  5 seconds,  which  gives  a = 32.2  ft/sec2 
and  k * 0.4  switch/sec. 

It  should  be  pointed  out  that  a system  design  based  on  the  target 
motions  described  above  would  operate  well  against  this  class  of  maneuver 
as  a whole.  This  appears  to  be  a desirable  procedure . Nevertheless , 
without  altering  the  parameters  this  system  would  not  be  expected  to 
operate  as  well  as  it  could  against  one  particular  target  maneuver  such 
as  a single  target  turn.  Even  then,  however,  it  can  be  shown  by  simu- 
lation studies  that  systems  optimized  for  the  statistical  maneuver  used 
herein  sure  essentially  optimum  for  the  single  turn  maneuver  as  well. 
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ANALYSIS,  RESULTS,  AND  DISCUSSION 
The  Wiener  Filter  Theory 


The  problem  of  minimizing  the  effects  of  noise  can  be  considered  to 
be  one  of  compromise.  At  one  extreme  for  which  the  system  response  is 
fast  the  error  becomes  large  because  of  the  ability  of  the  missile  to 
follow  the  noise  too  well.  At  the  other  extreme  for  which  the  system 
response  is  too  slow  the  error  becomes  excessive  because  of  the  difficulty 
the  missile  has  in  following  the  target  maneuvers . The  optimum  system  is 
one  which  will  compromise  these  two  situations  in  the  best  possible  man- 
ner. More  precisely  stated  the  problem  becomes:  Given  the  statistical 
characteristics  of  the  two  input  quantities,  target  motion  and  noise, 
what  is  the  optimum  transfer  function  which  will  minimize  the  mean-square 
error  The  answer  to  this  problem  can  be  determined  by  a statistical 

theory  known  as  the  Wiener  filter  theory.  The  final  result  of  this  theory 
is  an  integral  which  when  evaluated  represents  the  optimum  linear  transfer 
function.  This  transfer  franc t ion  is  given  by 


*o(i»>  - +7"T to  at  (3) 

2jt4+(w)  '-'o  j-oo  4 (a-) 

In  this  equation  4+  and  4"  have  the  following  meaning:  If  4 is  the 
spectral  density  defined  by  the  equation 


4 = 4m  + 4 


’N 


then  4+  and  4"  must  satisfy 

4+  4"  = 4 

where  4+  has  poles  and  zeroes  only  in  the  upper  half  of  the  complex  u> 
plane  and  4“  only  in  the  lower.  More  details  are  given  in  Appendix  A. 
The  derivation  of  equation  (3)  is  beyond  the  scope  of  this  report;  for 
this  derivation  see  references  4 and  10. 

The  transfer  function  Y0  of  equation  (3)  is  a mathematical  repre- 
sentation of  the  box  in  figure  1(b).  According  to  the  method  of  deriva- 
tion, the  transfer  function  must  be  physically  realizable  which  means 
that  the  system  is  not  required  to  respond  to  an  input  before  that  input 
occurs . 

It  might  be  pointed  out  that  a general  solution  exists  which  involves 
the  problem  of  prediction  as  well  as  filtering.  However,  since  it  is 
apparent  from  figure  1(a)  that  we  are  interested  only  in  the  missile 
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position  coinciding  with  the  present  target  position,  the  prediction  aspect 
has  "been  eliminated  from  the  general  solution  in  writing  equation  (3) . 

Restrictions  involved  in  the  theory.-  There  are  certain  restrictions 
implicit  in  the  derivation  of  equation  (3)  so  that  the  validity  of  its 
application  to  the  beam-rider  guidance  system  depends  on  how  well  these 
restrictions  are  met.  First,  the  input  quantities,  target  motion  and  noise, 
must  be  stationary  random  series  (see  ref.  1 for  a detailed  definition)  and 
be  defined  by  corresponding  power  spectra  which  are  continuous.  A discon- 
tinuity in  the  spectrum  might  be  due  to  a predictable  component  such  as  a 
sine  wave;  such  components  must  be  eliminated  from  the  input  before  the 
theory  can  be  applied.  It  is  generally  believed  that  displacements  at  the 
target  are  approximately  stationary  random  series  (ref.  3) . Since  the 
beam-rider  system  operates  from  these  displacements,  the  inputs  to  the 
beam- rider  system  are  also  stationary  random  series.  (In  contrast,  the 
inputs  to  a proportional-navigation  guidance  system  are  nonstationary  ran- 
dom series  since  angular  inputs  are  measured  by  the  missile  itself  and  the 
angles  tend  to  become  larger  as  the  range  decreases.)  Second,  the  solution 
is  based  on  linear  theory  and  furthermore  is  applicable  only  to  a system 
with  constant  coefficients  in  its  differential  equation.  On  the  other 
hand,  the  kinematic  loop  of  guidance  systems  generally  involves  a time- 
vari  ole  range  factor.  In  particular  for  the  beam-rider  system,  the  time- 
variable  factor  is  the  ratio  of  the  launcher-to-missile  and  the  launcher- 
to-target  ranges.  However,  since  the  miss  is  determined  primarily  by  what 
happens  near  the  end  of  flight,  during  which  the  variation  in  this  ratio 
is  small,  it  is  reasonable  to  assume  that  the  requirement  of  constant 
coefficients  is  approximately  met. 


Evaluation  of  the  optimum  transfer  function.-  Within  the  above 
restrictions,  it  is  possible  to  use  equation  (3)  to  evaluate  first,  the 
general  form  of  the  optimum  transfer  function  and  second,  the  numerical 
constants.  As  an  illustration  of  this  method  the  theory  has  been  applied 
with  certain  simplifying  assumptions  as  discussed  in  Appendix  A to  the 
target  motion  and  noise  characteristics  described  earlier.  For  this  case, 
then,  equation  (3)  can  be  evaluated  to  give  the  general  form  of  the  optimum 
closed- loop  and  corresponding  open- loop  transfer  functions  as  follows: 


(Tq2s2  + 2^s  + 1) 

(Tps  + 1)(T72s2  + 25rTrs  + 1) 


OO 


M0(s) 


= *4 


(l^s2  + 2kTas  + 1) 
s2(T^s  + 1) 


(5.) 


These  equations  can  be  evaluated  for  any  specific  case  and,  as  an  example , 
for  the  specific  target  motion  and  noise  characteristics  given  below 
equation  (2)  and  with  the  magnitude  of  the  noise,  N,  chosen  to  be 
15  ft2/radian/sec,  the  optimum  transfer  functions  become 
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(0.855  a2  + 1.41  8+1) 

(0.687  s + l) (0.490  s2  + 0.727  8 + 1) 


(6) 


h0(s) 


7 42  (0-855  s2  + 1.41  s + 1) 

'*  s2(2.50  s + 1) 


(7) 


The  chosen  value  of  N represents  a mid- value  between  the  expected 
extremes.  The  significance  of  this  choice  will  be  discussed  later. 

The  corresponding  frequency  characteristics  are  possibly  more  illus- 
trative to  the  control  designer.  The  solid  curves  in  figure  2(a)  show 
the  optimum  closed- loop  characteristic  for  the  above  conditions  and 
describe  the  characteristics  of  the  box  in  figure  l(b) . It  might  be  noted 
that  since  the  break  point  of  the  noise  spectrum  occurs  at  6 cps,  it  is 
essentially  flat  over  the  frequency  range  of  importance  of  the  optimum 
transfer  function.  Since  actual  systems  are  usually  designed  on  an  open- 
loop  basis,  the  corresponding  open- loop  characteristic  is  shown  in  fig- 
ure 2(b).  The  considerations  involved  in  achieving  the  characteristics 
of  figure  2 are  discussed  in  later  sections. 

RMS  error  performance.-  Although  the  Wiener  filter  theory  can  be  used 
to  define  the  optimum  transfer  function,  it  does  not  give  the  minimum 
error  directly.  This  must  be  evaluated  from  the  optimum  transfer  function. 
With  the  earlier  assumptions  as  to  the  form  of  the  target  motion  and  noise, 
reference  1 showB  that  the  total  mean- square  error  is  composed  of  target 
and  noise  components  which  can  be  evaluated  according  to  equation  (8) . 


I*  = 


f\CO 

Y0(iaj)  |2®ij(u>)du  + / |Y0(iw)  |24>u(u>)du> 


(8) 


The  quantity  Y0(iw)  is  obtained  from  equation  (4)  by  the  cubstitution 
s = iu.  It  is  possible  to  evaluate  these  integrals  mathematically  but 
expedience  has  shown  that  it  is  easier  and  more  instructive  to  use 
graphical  techniques.  For  the  example  under  discussion,  evaluation  of 
equation  (8)  gives: 

=13.7  feet 
n/ Cij2  = 6.7  feet 
Ve*  = 15.4  feet 
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The  later  figure  represents,  then,  the  theoretical  lower  limit  of  rms 
error  that  can  he  achieved  by  a linear  constant-coefficient  transfer 
function  corresponding  to  the  assumed  target  and  noise  characteristics. 

For  later  use,  the  integrands  of  equation  (8)  for  the  previous 
example  are  plotted  in  figure  3*  Physically,  these  curves  can  be  inter- 
preted as  the  frequency  distribution  of  the  energy  in  the  components  of 
error.  It  is  seen  that  the  error  energy  is  concentrated  in  the  frequency 
range  below  2 cps. 

As  mentioned  earlier,  it  is  important  to  consider  a large  range  of 
noise  magnitudes  for  a number  of  reasons.  For  example,  measurements  of 
glint  noise  are  subject  to  considerable  discrepancy  and  may  cover  a wide 
range  of  magnitudes.  Likewise,  the  noise  magnitude  may  vary  due  to 
different  sized  targets  or  different  attack  aspects.  In  order  to  assess 
these  effects  the  minimum  error  has  been  determined  as  a function  of  the 
zero  frequency  noise  magnitude  N of  equation  (l)  by  a procedure  identical 
to  that  above.  This  requires  the  determination  of  a different  transfer 
function  Y0  for  each  value  of  N,  as  in  Appendix  A for  N = 15  ft2/ 
radian/sec,  and  the  evaluation  of  the  resultant  error  by  equation  (8). 

The  result  is  shown  in  figure  4 by  curve  A.  Each  point  on  this  curve 
corresponds  to  a different  optimum  transfer  function  and  it  is  possible 
to  use  this  transfer  function  as  a guide  in  designing  the  guidance  system. 
The  value  of  this  curve  is  that  it  can  be  used  as  a standard  with  which 
to  compare  the  performance  of  any  system. 

Also  indicated  in  figure  4 is  the  operating  range  of  noise  magnitudes 
corresponding  to  the  class  of  targets  and  aspects  previously  mentioned. 

At  first  it  might  appear  that  each  noise  level  would  require  a change  in 
the  filter  or  system  so  as  to  maintain  optimum  performance.  However,  as 
indicated  by  curve  B,  if  a system  is  optimized  only  at  the  mid-value  of 
the  noise  range  shown  here  as  the  design  value,  the  performance  of  this 
system  for  other  noise  levels  will  deviate  to  a negligible  extent  from  the 
optimum  over  the  range  of  interest.  Thus,  to  obtain  near  optimum  perform- 
ance over  the  range  of  noise  magnitudes  likely  to  be  encountered  it  is 
necessary  to  optimize  only  for  the  design  value  of  noise.  This  is  a very 
fortunate  fact,  and  one  of  obvious  practical  importance.  In  a similar 
manner  it  can  be  shown  that  a change  in  the  break  point  of  the  noise 
spectrum  would  not  greatly  affect  the  minimum  error  curve  shown  in  figure  4 
as  long  as  the  noise  spectrum  is  essentially  flat  over  the  bandwidth  of  the 
optimum  transfer  function  of  figure  2. 

It  is  interesting  to  compare  these  results  with  the  performance 
obtained  by  disregarding  noise  theory  in  the  design,  that  is,  by  designing 
the  system  for  the  fastest  possible  response.  Two  examples  are  shown. 

Curve  C illustrates  the  performance  against  noise  which  is  obtained  for  a 
system  with  perfect  response  characteristics  (unity  transfer  function) . 
Curve  D illustrates  the  performance  against  noise  for  a more  realistic 
system  given  in  reference  11  (variable- incidence  missile).  For  this  case 
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the  important  saturating  elements  were  simulated,  and  the  system  was 
designed  for  the  fastest  possible  response  to  a step  in  the  beam  of 
100  feet.  The  comparison  with  the  optimum  of  the  curves  for  both  of  these 
example  systems  shows  that  the  performance  is  significantly  poorer  than 
the  optimum  performance  indicated  by  theory. 


Application  of  Wiener  Theory  to  the 
Beam-Rider  Guidance  System 


To  achieve  the  optimum  results  indicated  by  the  Wiener  theory  it  is 
necessary  to  design  the  guidance  system  with  frequency  characteristics 
approximating  the  optimum  specified  by  figure  2.  The  difficulties  in 
accomplishing  this  as  well  as  the  modifications  which  are  required  are 
discussed  in  the  following  sections. 


Limitations  in  the  application  of  the  Wiener  theory.-  There  are 
several  restrictions  on  the  possible  forms  of  the  transfer  function  which 
can  be  achieved  in  practice.  The  foremost  of  these  is  manifested  by  a 
consideration  of  the  accelerations  required  of  the  missile.  For  the 
optimum  system,  Y0,  the  mean  square  of  the  required  acceleration  is 

expressed  by 


r>°°  2 r00  2 

* / | Y0 ( ito)  | w4®T(u>)du > + / |Y0(iu>)  | u>44>N(u))dii) 


(9) 


It  is  easy  to  see  that  for  the  form  of  noise  spectrum  given  by  equation  (l) 
and  the  optimum  transfer  function  of  equation  (4)  the  spectral  density  of 
the  acceleration  required  by  the  noise  increases  with  frequency  at  the  high 
frequencies,  giving  rise  to  an  infinite  called-for  acceleration.  Obviously 
finite  values  of  acceleration  can  be  obtained  only  if  the  transfer  func- 
tion Y0  falls  off  at  high  frequencies  like  l/w2  or  some  greater  power. 

To  the  control  designer,  this  restriction  is  perhaps  more  readily 
interpreted  in  terms  of  a control-deflection  restriction.  Assuming,  for 
the  sake  of  the  argument , that  the  aerodynamic  transfer  function  &/yM 
were  everywhere  linear  (which,  of  course,  it  is  not,  since  practically, 
control  deflections  cannot  increase  without  limit) , the  analog  of  equa- 
tion (9)  would  become 


6* 


brjf 


poo  „ 2 p oo  1 2 

= / | Y0( iuj)  ] T^-(iw)  4T(u>)du>  + / |Y0(iw)|  ^-(itn)  | $N(w)dw  (10) 
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It  can  be  seen  from  equations  presented  in  later  sections  that  8/yw 
approaches  a constant  multiple  of  w2  at  high  frequencies.  Thus,  the 
conclusion  can  be  drawn  from  equation  (10)  that  finite  values  for  the 
control  deflection  will  be  obtained  for  any  practical  system  which  may  be 
designed  only  if  the  system  transfer  function  differs  from  the  optimum  at 
high  frequencies. 

At  first  it  would  appear  that  any  departure  of  the  system  transfer 
function  from  the  optimum  Y0  would  result  in  increased  rms  error  per- 
formance. However,  an  examination  of  the  error  spectrum  indicates  that, 
since  the  power  contained  in  frequencies  above  2 cps  is  small,  attenuating 
the  transfer  function  above  this  frequency  will  not  affect  the  error 
appreciably.  On  the  other  hand,  the  control  motion  spectrum  of  figure  ‘5 
shows  that  attenuation  in  the  transfer  function  above  2 cps  will  greatly 
reduce  the  required  control  motion  since  the  major  portion  of  control 
motion  power  is  concentrated  in  the  higher  frequencies.  From  these  con- 
siderations it  is  obviously  desirable  to  attenuate  the  transfer  function 
Y0  for  all  frequencies  above  the  bandwidth  of  appreciable  error  spectrum. 

The  attenuation  might  be  supplied  in  a variety  of  ways.  However, 
if  it  is  assumed  that  an  ideal  filter  could  be  inserted  in  the  system  so 
as  to  supply  infinite  attenuation  for  all  frequencies  above  2 cps  without 
affecting  the  response  to  frequencies  below  2 cps,  it  is  still  found  that 
for  missiles  similar  to  the  one  studied  here  deflections-  required  for  this 
frequency  band  are  larger  than  are  available.  Specifically,  for  the 
missile  used  herein  the  control  deflection  required  is  about  15°  rms; 
since  the  control  deflection  is  physically  limited  to  about  14°  (ref.  11), 
the  called- for  control  motions  are  larger  than  are  available.  Thus  it 
is  impossible  for  the  actual  control  motion  to  follow  the  required  pro- 
gram and  the  optimum  missile  motions  cannot  be  achieved.  Actually  the 
value  of  15°  is  quite  optimistic  since  the  ideal  filter  does  not  exist 
and  cannot  even  be  approached  very  closely.  Consequently,  the  simultane- 
ous achievement  of  the  optimum  error  spectrum  and  called-for  control 
motions  which  do  not  exceed  the  limited  values  are  conflicting  require- 
ments. The  importance  of  this  conflict  and  the  solution  of  the  problem 
posed  by  it  are  considered  in  the  following  sections. 

Optimum  filtering  with  limiting.-  Because  of  the  above  restrictions 
in  the  application  of  the  Wiener  theory,  it  is  generally  necessary  to  make 
certain  modifications  in  the  practical  case.  In  what  follows,  the  effects 
of  limiting  and  the  consequent  modifications  required  are  discussed. 

Effects  of  limiting:  There  are  many  limiting- type  nonlinearities 
which  may  exist  in  a guidance  system.  Those  of  the  most  importance  in 
the  beam-rider  system  are  (a)  control  motion  limiting  due  to  mechanical 
limitations,  (b)  rate  of  control  motion  limiting  due  to  restrictions  on 
servo  capability,  and  (c)  radar  receiver  voltage  limiting  due  to  cir- 
cuit restrictions.  Typical  values  which  are  used  in  these  studies  are 
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0.25  radian,  5 radians/sec,  and  100  volts,  respectively.  The  effects  of 
these  limits  can  he  severe  and  need  to  be  considered  in  the  filtering 
problem. 

By  means  of  Booton ' s recent  theory  (ref.  12)  it  is  possible  to  evalu- 
ate the  effects  of  nonlinearities  on  the  over-all  system  performance.  As 
applied  to  the  beam- rider  system,  the  theory  shows  that  it  is  possible  to 
approximate  each  of  the  saturable  elements  by  a simple  gain  which  can  be 
suitably  chosen  so  that  for  increased  rms  inputs  to  the  saturable  element 
this  gain  is  reduced.  Analysis  shows  that  in  terms  of  the  over-all  sys- 
tem characteristics  the  effects  of  this  gain  reduction  are:  (l)  to  reduce 
the  open- loop  system  gain  or  (2)  to  increase  the  time  constants  in  the 
open- loop  transfer  function.  In  general,  both  effects  result  in  a reduc- 
tion in  bandwidth  of  the  closed- loop  frequency  response.  In  this  respect, 
then,  limiting  effectively  results  in  additional  filtering  which,  of 
course,  extends  further  into  the  low-frequency  ranges  as  the  rms  inputs 
increase.  Thus  the  system  response  can  become  so  slow  as  to  result  in  a 
large  increase  in  error  due  to  the  missile's  inability  to  follow  the  tar- 
get motions.  Furthermore,  inspection  of  figure  2(b)  shows  that  the  Wiener 
transfer  function  is  conditionally  stable;  for  large  rms  inputs  the  Wiener 
transfer  function  may  become  unstable . That  this  does  occur  in  examples 
to  be  discussed  later  has  been  demonstrated  during  REAC  studies.  Conse- 
quently the  large  gain  reductions  and  phase  lags  introduced  by  limiting 
cannot  be  tolerated. 

Additior  j1  filtering:  In  general,  the  energy  of  the  saturating  quan- 
tities contained  in  the  higher  frequencies  is  much  larger  than  that  in  the 
lower  frequencies.  For  this  reason,  purposely  introducing  filtering  into 
the  system  so  as  to  attenuate  the  response  at  the  higher  frequencies  can 
be  effective  in  reducing  the  effects  of  limiting.  This  can  be  accomplished 
in  either  of  two  ways:  by  inserting  additional  networks  in  the  guidance 
system,  or  by  changing  those  inherent  time  constants  which  exist  in  any 
actual  system  and  which  are  not  required  for  the  optimum  filtering. 

Because  of  the  complexity  of  the  beam-rider  guidance  system  there  are  a 
great  number  of  places  where  the  additional  filtering  can  be  introduced. 
These  possibilities  are  discussed  later. 

The  selection  of  the  desired  frequency  characteristics  of  the  added 
filtering  is  more  difficult  than  the  selection  of  filtering  location.  In 
general,  the  choice  of  the  desired  amount  will  involve  a compromise  between 
two  extremes:  (l)  The  first  extreme  is  one  in  which  the  added  filtering  is 
too  severe.  As  the  added  filtering  is  extended  into  the  lower  frequencies, 
limiting  effects  are  reduced  but,  as  previously  mentioned,  the  results  of 
limiting  can  only  be  eliminated  at  the  expense  of  increased  error.  (2)  The 
other  extreme  is  one  in  which  the  added  filtering  is  insufficient  to  avoid 
serious  limiting  effects.  This  extreme  is  similar  to  the  previous  one 
inasmuch  as  the  effects  of  limiting  are  also  equivalent  to  additional  fil- 
tering. For  serious  limiting,  this  equivalent  filtering  is  too  severe  and 
results  in  increased  error  and  possibly  instability. 
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It  is  apparent  that  the  effects  of  .Limiting  and  of  added  filtering 
by  means  of  networks  are  closely  related,  in  that  both  are  equivalent  to 
the  introduction  of  additional  filtering  into  the  system.  Extremes  of 
either  type  reduce  the  performance  of  the  system  sufficiently  that  the 
error  increases  due  to  the  inability  of  the  missile  to  follow  the  target. 


Design  of  the  Guidance  System 


In  applying  the  modifications  discussed  in  the  preceding  paragraphs 
to  the  design  of  the  guidance  system  it  is  found  convenient  to  consider 
two  separate  stages:  The  first  stage  involves  the  design  of  the  system 
on  a linear  basis  to  approximate  the  optimum  Wiener  transfer  function, 
while  the  second  involves  supplying  the  additional  filtering  terms  to 
minimize  the  limiting  effects.  The  optimum  Wiener  transfer  function  may 
be  designed  into  the  missile- control  system,  or  into  the  tracking  radar, 
or  may  Le  apportioned  between  the  two,  with  the  additional  filtering  in 
either  place.  Thus,  there  are  a number  of  possibilities.  Obviously,  for 
linear  systems  all  approaches  could  be  designed  so  as  to  produce  identical 
results.  In  general,  this  does  not  hold  for  nonlinear  systems,  that  is, 
systems  in  which  limiting  occurs,  so  that  in  the  practical  case  it  is 
necessary  to  consider  the  effect  of  the  filtering  location  on  performance. 
The  relative  merits  of  placing  the  optimum  Wiener  filtering  first  in  the 
missile-control  system  and  second  in  the  tracking  radar  are  discussed  in 
the  next  two  sections. 

Missile-control  system  designed  for  optimum  Wiener  filtering.-  For 
the  case  of  missile-control-system  filtering  it  is  necessary  that  the 
tracking  radar  be  designed  fast  enough  to  follow  the  target  motion  and 
the  noise,  since  the  missile- control  system  performs  the  optimum  filter- 
ing operation.  The  first  step  in  the  design  of  the  missile-control  system 
is  the  synthesis  of  a linear  system  which  matches  the  optimum  Wiener 
transfer  function.  It  has  been  shown  that  it  is  only  necessary  to  approx- 
imate this  transfer  function  over  the  frequency  range  of  appreciable  error. 
This  can  be  accomplished  in  many  ways,  but  to  achieve  this  design  by  a 
system  of  conventional  form,  the  design  is  most  easily  by  cut- and- try 
procedures . In  this  procedure  it  is  desirable  to  select  the  aerodynamics 
as  a starting  point  since  the  design  of  the  airframe  is  relatively  inflexi- 
ble compared  to  the  design  of  the  control  system.  The  characteristics  of 
a typical  variable-incidence  missile  were  chosen  from  reference  11.  The 
transfer  functions  in  the  yaw  plane  for  this  missile  are  given  by  the  fol- 
lowing equations: 


+ 1 

Tds(Ta2s2  + 2?aTas  + 1) 
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lU  Vs2  4-  2ghTbS  -l-  1 

5 * Ts2s2(Ta2s2  + 2£aTas  + 1) 

Table  I summarizes  the  values  of  the  parameters  for  this  missile  for  a 
given  operating  condition  which  is  held  fixed  in  this  report.-  For  this 
missile  it  is  possible  to  synthesize  a linear  system  to  approximate 
closely  the  optimum  transfer  function.  One  possible  system  is  illustrated 
in  figure  6 and  the  corresponding  over-all  open- loop  transfer  function  is 
given  by 


kgkpV  (Tips  4 1)(T4s  + l)(Tb2s2  + g^TfaS  + 1) 

^ Td  + xsk3  s2(TnS  + l)(T5s  + l)(a83  + bs2  + cs  + l)  1 

The  derivation  of  this  transfer  function  is  given  in  Appendix  B.  It 
should  be  noted  that  the  transfer  function  of  equation  (ll)  differs  from 
that  of  equation  (5) ; however,  by  choosing  the  parameters  as  given  in 
Appendix  B and  as  tabulated  in  table  II,  column  G>,  the  two  transfer 
functions  can  be  closely  matched  over  the  frequency  range  of  interest  as 
illustrated  by  the  dotted  curves  of  figures  2 and  3*  It  might  be  noted 
that  the  only  restriction  on  the  parameter  T5  in  order  to  keep  the  two 
transfer  functions  closely  matched  is  that  it  be  small.  Thus  T5  may 
be  varied  somewhat  without  much  penalty. 


The  particular  design  discussed  above  has  been  chosen  only  for  optimum 
noise  performance  so  that  other  considerations  important  in  an  over-all 
design  might  dictate  certain  modifications.  These  possibilities  are  dis- 
cussed later.  It  is  also  apparent  that  since  the  design  of  the  above 
system  has  been  based  on  assumptions  of  linearity  the  performance  indicated 
above  cannot  be  achieved  in  practice  due  to  limiting  effects.  As  indicated 
earlier  these  limiting  effects  may  result  in  instability.  Hence  modifica- 
tions are  required.  These  modifications  in  the  two  cases  of  additional 
filtering,  first  in  the  tracking  radar  and  second  in  the  missile-control 
system,  are  discussed  in  the  following  paragraphs. 

Additional  filtering  in  tracking  radar:  The  simpler  means  ^f  intro- 
ducing additional  filtering  into  the  guidance  system  in  order  to  reduce 
limiting  is  to  place  it  in  the  tracking  radar  since  the  added  filtering 
and  optimum  Wiener  filtering  are  achieved  separately  by  the  tracking  radar 
and  missile-control  system,  respectively.  This  separation  has  the  advan- 
tage of  allowing  the  additional  filtering  to  be  altered  without  affecting 
the  optimum  filter  design  of  the  missile -control  system. 

The  additional  filtering  in  the  tracking  radar  may  take  innumerable 
forms.  For  this  study  the  closed- loop  transfer  function  of  the  tracking 
radar  was  assumed  to  be  of  the  following  form: 

^B  TiS  + 1 

= fxg«g  + TXs  4 1 (12) 
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This  transfer  function  was  chosen  to  provide  both  simplicity  and  a zero- 
velocity  error  system  and  it  can  be  approximated  by  properly  shaping  the 
networks  in  the  tracking  radar.  With  this  system,  the  amount  of  filtering 
can  conveniently  be  varied  by  a change  in  the  time  constant  ti.  At  the 
one  extreme  for  which  ti  is  very  small  the  filtering  provided  by  the 
tracking  radar  is  negligible,  thus  resulting  in  an  unstable  missile-control 
system  because  of  limiting  effects.  As  the  time  constant  of  the  radar  is 
increased,  however,  more  filtering  is  provided  and  limiting  effects  are 
reduced.  By  varying  figure  7 was  obtained;  it  shows  the  variation 
of  rms  error  with  the  natural  frequency,  1/2«ti,  of  the  tracking  radar. 

For  comparison,  the  minimum  error  obtained  from  Wiener  theory  is  also 
shown.  As  would  be  expected  from  previous  discussion,  the  optimum  opera- 
ting point  occurs  at  a tracking  radar  frequency  which  is  within  the  pass- 
band  of  the  missile-control  system.  At  this  frequency  limiting  effects 
still  exist  but  are  not  too  serious.  The  increase  in  error  above  the 
Wiener  theory  result  can  be  attributed  to  the  additional  filtering  intr<. 
duced  in  the  tracking  radar  and  to  the  limiting  effects  in  the  missile 
control  system.  The  rms  error  performance  of  the  system  defined  by  the 
optimum  operating  point  is  given  by  curve  E in  figure  8 against  a variable 
noise  level. 

Additional  filtering  in  missile-control  system;  An  alternative  place 
to  introduce  additional  filtering  into  the  guidance  is  the  missile-control 
system.  This  could  be  introduced  by  an  additional  network  in  the  radar 
receiver.  However,  since  it  has  been  pointed  out  that  the  choice  of  time 
constant  T5  was  somewhat  arbitrary,  it  is  possible  to  increase  T5  to 
provide  this  additional  filtering.  Thus  limiting  effects  can  be  reduced. 

The  introduction  of  added  filtering  into  the  missile-control  system 
has  the  disadvantage  that  the  added  filtering  also  affects  the  system 
stability  since  an  increase  in  a time  lag  is  destabilizing.  Hence,  added 
stability  is  required.  This  could  be  provided  by  altering  any  of  the 
basic  parameters  which  are  responsible  for  the  conditional  stability 
characteristic  of  the  Wiener  transfer  function.  To  preserve  the  low- 
frequency  characteristics  of  the  optimum  system  as  few  changes  as  possible 
are  desired.  The  parameter  Tn  is  a convenient  one  with  which  to  intro- 
duce this  stability. 

The  effects  of  the  above  changes  were  investigated  by  means  of  analog 
simulation.  Since  in  the  case  of  missile-control-system  filtering  the 
tracking  radar  is  not  required  to  filter,  it  should  be  designed  to  respond 
quickly.  For  this  purpose  the  radar  transfer  function  of  the  form  of 
equation  (12)  was  utilized  by  choosing  the  constants  so  as  not  appreciably 
to  alter  the  input  spectra.  This  was  accomplished  by  the  choice  of  a 
tracking  radar  natural  frequency  of  6 cps.  The  parameters  in  the  missile- 
control  system  were  altered  according' to  the  above  discussion.  Typical 
results  obtained  for  the  system  optimized  only  at  the  design  value  of 
noise  are  illustrated  by  curve  F in  figure  8.  The  optimum  Wiener  perform- 
ance is  repeated  here.  In  general,  results  similar  to  curve  F can  be 
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obtained  by  several  different  combinations  of  parameters . One  of  these 
combinations  requiring  the  fewest  modifications  from  the  Wiener  system  is 
given  by  the  parameters  in  column  © of  table  II.  It  can  be  seen  that 
the  additional  filtering  is  supplied  by  an  increased  radar  receiver  net- 
work time  constant  T5  while  the  added  stability  is  introduced  by  a 
decrease  in  the  time  constant  Tn.  In  this  case  the  increase  in  error 
above  the  optimum  is  due  to  the  added  filtering  introduced  by  the  network 
and  limiting,  and  to  the  slight  alteration  in  the  Wiener  system  time  con- 
stant T11. 

Tracking  radar  designed  for  optimum  Wiener  filtering.-  In  this  section 
are  considered  the  results  obtainable  when  all  of  the  optimum  filtering  is 
located  in  the  tracking  radar.  Here  the  desired  optimum  tracking  radar 
transfer  functions  are  given  by  equations  (4)  and  (5),  and  ideally  the 
missile-control  system  should  have  a transfer  function  of  unity.  However, 
for  many  reasons,  principally  those  arising  from  nonlinear  effects,  the 
latter  may  be  expected  to  depart  considerably  from  the  ideal.  Consequently 
it  is  desirable  to  design  the  missile-control  system  for  the  fastest  pos- 
sible transient  response  within  the  limitations  of  these  nonlinearities. 

Two  such  missile-control  systems  were  considered.  The  first  was  a con- 
ventional system  chosen  from  reference  11  in  which  the  response  time  to 
a step  of  100  feet  in  the  beam  was  minimized.  The  second  was  a system 
optimized  for  minimum  response  time  to  a small  enough  step  so  that  lin- 
earity was  not  exceeded.  The  control- system  parameters  for  these  systems 
are  given  in  columns  0 and  0 , respectively,  in  table  II. 

The  rms  error  performance  obtainable  for  tracking -radar  filtering 
is  summarized  in  figure  8 by  curve  G.  It  was  found  that  both  the 
missile-control  systems  gave  essentially  the  same  results.  In  this  case 
the  increase  in  error  above  the  Wiener  optimum  is  due  to  slowness  of  the 
missile-control  system  and  its  failure  to  follow  the  beam  perfectly 
because  of  limiting  effects. 

It  might  appear  that  these  limiting  effects  could  be  reduced  by  addi- 
tional filtering  in  the  tracking  radar.  This  possibility  was  explored  by 
the  addition  of  a simple  first-order  filter  to  the  optimum  open-loop  trans- 
fer function  given  in  equation  (5).  By  increasing  the  time  constant  of 
the  added  filter  the  tracking -radar  performance  is  deteriorated  from  the 
theoretical  optimum  while  the  error  of  the  missile-control  system  is 
decreased  because  of  a reduction  of  limiting  effects.  Results  of  these 
combined  effects  are  shown  in  figure  9 from  which  it  can  be  seen  that  the 
added  filtering  results  in  progressively  poorer  over-all  performance.  Thus 
the  beneficial  effects  of  decreased  limiting  are  overbalanced  by  the  detri- 
mental effects  of  altering  the  tracking -radar  transfer  function  from  the 
optimum. 
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Comparisons  and  Other  Considerations 


The  results  of  various  filtering  arrangements  are  compared  in  fig- 
ure 8 on  the  basis  of  the  rms  error.  If  it  were  not  for  the  nonlinearities 
the  performance  for  both  missile-control  system  filtering  and  tracking- 
radar  filtering  could  be  made  identical  with  the  performance  of  the  optimum 
Wiener  system.  However,  it  is  interesting  to  find  that  even  in  the  non- 
linear case  all  arrangements  can  result  in  comparable  error  performance, 
dome  advantage  in  tracking- radar  filtering  is  apparent  from  the  figure. 

For  comparison,  again,  the  performance  of  a typical  missile-control  system 
which  was  optimized  for  fast  response  in  the  absence  of  noise  is  repeated. 

Up  to  this  point  only  the  error  performance  has  been  considered. 

Other  factors  which  are  of  importance  in  the  over-all  evaluation  of  a guid- 
ance system  will  now  be  discussed.  One  such  factor  is  the  servo  energy 
required  to  achieve  a given  error  performance.  For  a servo  system  already 
designed  to  meet  the  maximum  expected  hinge  moment,  the  servo  power  is 
proportional  to  the  time  average  of  the  sum  of  the  absolute  displacements 
of  S between  values  at  which  5 changes  sign.  Thus  the  average  servo 
power  for  a time  t is 


ZlA&l 

average  servo  power  ~ — ^ — L 

Evaluation  of  the  servo  power  for  the  systems  discussed  above  has  shown 
that  both  methods  of  missile-control-system  filtering  require  about 
22  percent  more  power  than  does  tracking- radar  filtering.  Thus  in  an 
over-all  evaluation  based  both  on  rms  error  performance  and  servo  energy 
requirements  it  is  apparent  that  tracking  radar  filtering  is  slightly 
superior . 

There  may  be  still  other  factors  of  importance  in  guidance  system 
design  even  within  the  framework  of  the  assumptions  discussed  earlier. 
Usually  these  requirements  are  related  to  the  specific  design  objectives 
of  the  system  and  may  dictate  certain  modifications  such  as  the  choice  of 
filtering  location  or  alterations  in  certain  individual  transfer  functions 
of  the  system.  For  example,  for  short-range  missiles  in  which  launching 
errors  are  not  prevented  from  building  up,  capture  of  the  beam  in  minimum 
time  may  be  important  enough  that  tracking-i odar  filtering  would  be  pref- 
erable to  missile-control-system  filtering  because  of  the  fast  response 
which  can  be  designed  into  the  missile-control  system.  In  other  cases 
it  may  be  desirable  to  alter  the  design  of  the  system  somewhat  for  any  of 
a number  of  reasons.  Requirements  of  simplicity  on  certain  parts  of  the 
system  or  the  necessity  of  using  certain  fixed  and  unalterable  elements  in 
rhe  design  are  examples.  Another  possibility  is  that  for  flight  conditions 
in  which  serious  atmospheric  turbulence  exists  it  might  be  necessary  to 
make  certain  alterations  to  minimize  the  response  to  gusts.  An  investiga- 
tion of  such  factors  is  beyond  the  scope  of  this  report. 
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CONCLUDING  REMARKS 

This  study  has  considered  the  application  of  Wiener  filter  theory  to 
the  optimization  of  a beam-rider  system  operating  in  the  presence  of  noise. 
The  Wiener  theory  has  been  used  to  establish  both  the  theoretical  lower 
limit  of  error  and  the  desired  transfer-function  characteristics.  Although 
these  transfer-function  characteristics  might  be  expected  to  vary  with 
noise  magnitude,  it  has  been  found  that  when  the  transfer  function  is  cho- 
sen to  be  optimum  for  only  a mid-value  of  noise  the  performance  for  other 
noise  levels  deviates  to  a negligible  extent  from  the  optimum. 

In  general,  linear  beam-rider  systems  can  be  synthesized  to  produce 
theoretical  optimum  performance,  but  the  practical  achievement  of  these 
results  is  restricted  by  limiting  of  control  deflections,  rate  of  control 
deflections,  and  radar  receiver  voltages.  With  suitable  modifications, 
however,  the  theory  can  be  useful  as  a guide  in  system  design.  Optimiza- 
tion in  the  actual  nonlinear  ca3e  is  shown  to  involve  two  considerations: 
limiting  and  additional  filtering.  It  is  shown  that  the  important  types 
of  limiting  tend  to  result  in  system  instability.  However  by  appropriate 
placement  of  additional  network  filtering  it  is  possible  to  minimize  these 
limiting  effects  without  serious  deterioration  of  the  error  performance. 

The  design  of  a guidance  system  is  most  conveniently  accomplished  in 
two  stages:  The  first  consists  of  designing  the  system  on  a linear  basis 
to  approximate  the  Wiener  transfer  function,  and  the  second  of  supplying 
the  additional  filtering  terms  to  minimize  the  limiting  effects.  The 
application  of  the  Wiener  theory  and  the  modifications  required  to  arrive 
at  an  optimum  system  design  have  been  illustrated  in  this  report  by  con- 
sidering systems  in  which  the  optimum  Wiener  filtering  is  designed  into 
either  the  missile-control  system  or  the  tracking  radar,  and  additional 
network  filtering  supplied  in  either  place.  It  was  found  that  comparable 
error  performance  can  be  achieved  by  any  of  the  methods.  Consideration 
of  both  rms  error  and  servo  energy  requirements  for  the  cases  studied 
indicates  that  optimum  filtering  in  the  tracking  radar  is  slightly  superior 
to  that  in  the  missile-control  system. 

In  modifying  the  results  of  the  Wiener  theory  in  order  to  minimize 
the  effects  of  limiting,  the  best  results  were,  in  general,  obtained  when 
the  additional  filtering  was  added  in  such  a way  as  to  tend  to  keep  the 
operation  of  the  system  in  the  linear  range.  This  has  suggested  the  pos- 
sibility of  seeking  an  optimum  solution  based  on  the  stipulation  that  the 
filtering  should  restrict  the  operation  of  the  system  to  within  its  linear 
range.  In  reference  13  the  results  of  such  an  analysis  are  presented  and 
compared  to  those  of  the  present  report. 
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APPENDIX  A 


DETERMINATION  OF  THE  OPTIMUM  WIENER  FILTER 


Wiener’s  solution  for  the  optimum  transfer  function 
expressed  by  the  following  equation: 

Tod., . -krt  r*-iut  r ** « 

2jt®+(ai)  J0  JLoo  $ (a) 


where  represents  the  target  motion  spectral  density,  and  the  quanti- 

ties ®+  and  9~  are  defined  as  the  factors  of  a certain  function  <b  with 
poles  and  zeroes  located  in  the  upper  and  lower  half -planes,  respectively; 
the  factor  <f  is  defined  as  equal  to  so  that 

9 = ®rp  + (A2) 

For  the  case  to  be  considered  here,  the  target  motion  and  noise  are 
defined  by 

k2)  (A3) 

»H  = H (A4) 

It  will  be  noted  that  the  noise  spectrum  has  been  approximated  here  by  a 
constant  in  order  to  reduce  the  complexity  of  the  calculations.  This 
approximation  is  valid  because  the  more  exact  form  of  the  noise  spectrum 
(eq.  (l))  is  essentially  flat  over  the  bandwidth  of  the  optimum  transfer 
function.  Use  of  the  more  exact  form  of  the  noise  spectrum  affects  only 
the  response  at  the  high  frequencies  which  are  beyond  the  range  of  inter- 
est. It  should  also  be  pointed  out  that  the  use  of  equation  (A3),  as 
such,  leads  to  certain  mathematical  difficulties  in  evaluating  the  right- 
hand  side  of  equation  (Al),  because  the  theory  requires  that  the  poles 
of  not  be  located  on  the  real  axis.  Rigorously,  to  avoid  these 
diff icurties , it  is  necessary  to  modify  the  target  motion  spectral  density 
to  the  following: 


— 


rt(u)2  + Hi2)^  + + k2) 


where  rji  and  rj2  are  any  small  real  numbers.  Thus  the  solution  of  equa- 
tion (Al)  is  a function  of  iw,  rji,  and  rjg.  The  desired  answer  is  then 
obtained  as  the  limit  of  Yo(iw,Tli/Tte)  612  Hi  and  T|2  approach  zero. 
However,  it  can  be  showr  that  the  same  answer  can  be  obtained  more  simply 
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by  taking  t]i  and  r\2  equal  to  zero  as  would  be  obtained  fran  the  rigor- 
ous process  described  above.  It  will  be  convenient  to  retain  the  tj’s 
for  a few  lines;  subsequently  they  will  be  dropped.  According  to  equa- 
tion (A2)  then, 


. . ae  + k^a4  + (ka2/itN) 

a T + (a  + irji)(a  + irj2)(a  - irji)(a  - iTig)(a2  + k2) 

(A6) 

a>(„\  - w (a  ~ ai)  (a  - ag)  (a  - dp)  (a  + ai)  (a  + ag)  (a  + op)  , , 

(a  + itji)(a  + irj2)(a  - irji)(a  - itj2)(a  + ik)  (a  - ik) 


where  ax,  az,  and  03  represent  roots  in  the  upper  half -plane.  It  should 
be  noted  that  none  of  these  roots  are  real.  From  equation  (A2)  also,  it 
follows  that 


$+(a)  = » (a  ~ ai)(a  ■ ge)(a  - oq)  _ N p(a) 

(a  - it)i)(a  - iiteHa,  - ik)  (a  - Itjx)  (a  - itjgHa  - ik) 


= (a  + ai)(a  + Og)(a  + 03)  ^ -p(-q)  

' ’ (a  + irji)  (a  + 1^2)  (a  + ik)  (a  + it)i)(a  + iT)2)(a.  + ik) 


fA8) 


(A9) 


According  to  the  definition  given  above 

p(a)  = (a  - ai)  (a  - ag)  (a  - 03) 
3 a3  + bga2  + bia  + b0 


(A10) 


and 


p(-a)  = a3  - b^2  + bia  - bQ 


where 

b2  = -(ax  + a®  + 03 ) 
b x = aide  + aids  + ogae 

b©  * -axa^ao 


(All) 


(A12) 
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The  inner  integral  of  equation  (Al)  is 

Mtt)elatda  _ ka2  f°° && ^ 

4>"(<x)  wi*  a2  (a  - ik)(a  + ai)(a  + cte)(a  + 03) 

= ^ir  j^  f(a)dja  (A13) 

In  the  above  expression  it  should  be  noted  that  it  has  been  assumed  t)i 
and  tj2  are  zero  in  order  to  simplify  the  following  expressions.  Now  if 
a is  considered  to  be  a complex  variable,  the  above  integration  is 
equivalent  to  integration  over  the  contour  shown  below. 


Here  the  only  two  poles  involved  are  a second-order  pole  at  the  origin 
(actually  at  irji  and  irjg  in  the  rigorous  case)  and  a first-order  pole 
at  a = ik  since,  as  indicated  above,  -a-i,  -02,  and  -03  lie  in  the 
lower  half -plane . The  pertinent  residues  can  be  found  in  the  usual  man- 
ner as  shown  below: 


Res(o)  = lim  [a2f(a)] 
a o 


* + JL_  - b2 

kb0  k^o2 


(A14) 


Res(ik)  = lim  [(a  - ik)f(a)] 
a ■>  ik 


e~kt e-kt 

-1c2 (ik  + ai)(ik  + 02)  (ik  + <33)  -k2[-p(-ik)] 


(A15) 
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Hence  the  integral  in  equation  (A13)  becomes 


Ii 


= 2ka2ii  - 


= 2ka2i 


7i 


where  the  definitions  of  y 


ikbi  + bo  t e-kt  1 

kSbo2  kbo  -k2[-p(-ik) ]J 

+ Efe  + ^e'kt] 

(Al6) 

x and  y2  are  apparent. 

The  second  integration  in  equation  (Al)  is  merely  a Fourier  transform 
and  is  found  as  follows: 


XOO 

e-iwtiidt 


/00i- 

yi6-iu)t  + e-ia)t  + + iw)t  (A17) 


Thus, 


l2  = 2ka2i  + 


\h 

|_iu> 


kb0(iu)2  + k 


i-1 

+ iwj 


= 2ka2i 


(7l  + 7a)(iw)2  + (k  7l  + flQim  + ^ 
(iw)2(k  + iu>) 


(A18) 


It  should  be  noted  that  because  Tji  and  r\2  were  assumed  to  be  zero  for 
simplicity,  questions  concerning  the  existence  of  the  above  integral 
arise.  However  if  the  analysis  is  made  without  this  assumption,  it  can 
be  shown  that  the  integral  in  equation  (A17)  does  exist,  and  that  the 
limit  of  this  integral  as  t)i  ^ H2  so  bo  zero  becomes  precisely  equa- 
tion (Al8) . The  coefficients  in  this  equation  can  be  simplified  to  a 
more  useful  form  by  the  following  development.  From  the  definitions 
given  earlier  in  equations  (A6)  through  (A9), 

a6  + k2^  = [p(a)][-p(-a)] 


which  at  a = ik  becomes 


-gj'  = -b02  = [p(ik)  ] [-p(-ik)  ] 
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When  this  relation  is  substituted  in  the  definition  of  given  in 
equation  (Al6) , 


7 2 


-k?[  -p(-ik)  ] 


P.t1*) 

(-k2)(-b02) 


-ik3  - bak2  + ibjk  + b0 
^ 


Thus  the  coefficients  in  equation  (Al8)  simplify  to 


7i  + 7Z 


ik  + b2 


tyi  + 


The  transfer  function  Yc 
equations  as 


can  now  be  found  from  the  preceding 


Y0(iu)) 


la 

2 «®+ 


- **&*  (*>•  - y w * 1 
(-  # + & + 0(-  & + 0 


(A19) 


which,  in  terms  of  the  conventional  complex  frequency  s = iu>,  reduces  to 
the  following  alternative  forms: 


Y0(s) 


(#  • * >X#  * *) 


Y0(s) 


Y0(s)  = 


(-  Hr*)* : (•  & - , 

Tq2s2  + 2^Tns  + 1 


(Tps  + 1)(T72s2  + 2£7T7s  + 1) 


(A20) 
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where 


rr  2 _ . 

ia  - 


ik  + bg 
bo 


o&oq 


2?aTa 


ibi 


153 


2hT7 


. ag  + aa 
CtcOQ 


In  terms  of  the  equivalent  open-loop  transfer  function,  n0,  of  a unity- 
feedback  system, 


where 


Ho(8) 


Yq(s) 

1 - Yo(fl) 


ki 

ik 


(- 


+ 1 


* 1^ 


Tq282  + 2£aTqS  4-  1 
82(TxS  + 1) 


(A21) 


The  coefficients  of  the  transfer  function  of  the  optimum  system  have 
been  evaluated  for  the  following  values  of  target  motion  and  noise: 

k = 0.4 

$ = 5 sec 
a = 1 g 

N = 15  ft2/radian/sec 

Evaluation  of  the  numerator  of  equation  (A6)  gives 
a®  + k^a4  + = a®  + 0.l6  a4  + 8.8 

= (a  - ai)(a  - az)(a  - 03)  (a  + a,i)(a  + aa)(a  + 03) 
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where 


cl\  = i 1.45  6 


02  = 1.427  exp(i  0.543) 
oq  = -1.427  exp(-i  0.543) 

which  are  all  located  in  the  upper  half -plane.  These  roots  are  then 
sufficient  to  determine  the  optimum  transfer  function.  The  constants  in 
equation  (A12)  give 


h2  = -i  2.94 
bi  = -4.201 
b0  = i 2.97 

Then  evaluating  the  parameters  below  equation  (A20)  and  in  equation  (A21) 
gives  the  following  values: 


Parameter 

Value 

13a 

0.925 

5a 

.765 

t3 

.687 

t7 

.700 

£7 

.519 

Tx 

2.50 

7.42 

Fran  these  parameters  the  optimum  transfer  function  of  equation  (A20) 
becomes 


Y (b) (0.855  b2  + 1.41  s + 1) 

0V  ' " (0.687  S + 1) (0.490  s2  + 0.727  S + 1) 


(A22) 


and  from  equation  (A2l)  the  equivalent  open- loop  transfer  function  becomes 


ho(s) 


7 kp  (0.855  a2  + 1.41  s + 1) 

s2(2.50  s + 1) 


(A23) 
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APPENDIX  B 

MISSILE-CONTROL- SYSTEM  APPROXIMATION  TO  THE  OPTIMUM 
TRANSFER  FUNCTION 


An  optimum  linear  design  is  illustrated  in  figure  6,  and  in  the  fol- 
lowing the  system  equations  are  derived.  From  figure  6, 


It  (s>  = 

% 


ks(l^+l)(Tas+l) 


TdsCTxs+l)  (Ta2s2+2CaTa8+l)  (Tas+D+kskasC^js+l)  (T2s+l) 


(Bl) 


To  simplify  equation  (Bl)  it  is  convenient  to  make  T3  = T^.  This  choice 
is  not  essential  but  its  use  leads  to  more  easily  handled  equations.  In 
certain  cases  where  gust  disturbances  are  serious  it  may  be  more  desirable 
to  choose  T3  small.  With  the  former  assumption,  however. 


4<b) 


kB(TMS+l) 


s[  ( TdTa2Tx ) s3+  ( TdTa2+Td2 £aTaTi ) S2+  (TdT1+Td2 £aTa+ k8k3T2 ) s+  ( Td+k8k3 ) ] 


h§_ 


•%s+l 


Td+ksk3  s (as3+b82+c8+l) 
where 


a = 


b = 


TdTa2Tx 
Td  + ksk3 

Td(Ta2  + 2$aTaTx) 
Td  + ksk3 


_ ggaTaTd  + TxTd  + ksk3T2 
C Td  + kgka 

* 

Thus  the  entire  open-loop  transfer  function  can  be  written  as 

yM  kgk8v  (T4S  + 1) (Tips  + 1) (Tb2s2  + 2^bTbs  + l) 

^ 5 i'd  + ^8^3  s2(T5s  + l)(TnS  + l)(as3  + bs2  + cs  4 1) 


(B2) 


(B3) 


W 


It  is  desirable  to  choose  the  parameters  in  this  equation  so  as  to  match 
the  optimum  transfer  function  given  in  equation  (A2l) . Since  the  system 
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equation  (B4)  is  of  much  higher  order  than  the  optimum,  it  is  apparent 
that  the  matching  cannot  be  accomplished  perfectly.  However,  as  shown 
in  the  text  it  is  only  necessary  to  approximate  the  optimum  transfer  func- 
tion over  a limited  frequency  band.  This  can  be  accomplished  by  choosing 
certain  terms  in  (B4)  to  correspond  to  terms  in  (A21).  The  remaining 
terms  should  then  be  chosen  to  have  negligible  effect  on  the  optimum  error 
spectrum. 

There  are  many  ways  the  design  can  be  accomplished.  No  attempt  will 
be  made  to  investigate  all  the  possibilities.  One  design,  however,  is 
based  on  the  following  correspondences: 


T4  <-»  Ta 
Tio  <■'—>  Ta 
Tn  > Tx 


The  remaining  terms  can  be  chosen  in  any  manner  as  long  as  their  effect 
is  small  over  the  optimum  error  spectrum.  One  of  the  many  possible 
choices  is  to  factor  the  cubic  to  approximate  the  aerodynamic  term  in  the 
numerator  as  follows: 


as3  + bs2  + C8  + 1 = (Tbas2  + 2^-jrjTbs  + l)(TLs  + l)  (B5) 


where  TL  is  arbitrarily  chosen  to  be  small . This  factoring  is  not 
essential  to  the  design  but  it  leads  to  simple  equations  for  the  control- 
system  parameters.  For  example,  for  given  aerodynamics  and  choice  of 
Tl  the  coefficients  of  the  left  side  of  (B5)  are  determined  by 


a - TLTb2 
b = TL2^Tb  + Tb2 
c = tl  + 25bTb 


(B 6) 


Solving  equation  (B3)  then  for  the  control- system  parameters  gives 


Tl 

2?bTl,TLTa  + Tba(Ta  - 2CaTL)  L 

- Td  ■ 

2Th2 

Ta  ■ Tad  - ' 5aT»)  + T*  + 


> 


J 


(B7) 
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Then  the  value  of  kgk8  is  determined  from  the  desired  gain: 


k2k8 


Td  + kska 
— — v 


It  will  he  noted  that  since  only  k2k8  and  k8ks  are  specified,  one  of 
the  three  may  be  chosen  arbitrarily.  The  numerical  values  of  the  above 
parameters  have  been  calculated  and  are  listed  in  column  (^)  , table  II. 
It  should  be  pointed  out  that  this  design  is  only  one  of  many  possible 
designs.  In  any  particular  case  certain  modifications  may  be  necessary 
as  discussed  in  the  section  on  comparisons  and  other  considerations. 
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Figure  2.-  Optimum  transfer  functions. 
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(b)  Open-loop  transfer  function, 


Figure  2.-  Concluded. 
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Figure  3.-  Spectral  density  of  error  components. 
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Noise  spectral  density  magnitude,  N,  ft2/ radian /sec 

Figure  4.-  Optimum  performance  by  Wiener  theory. 
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Missile- control  system 


Figure  6.-  Linear-system  approximation  of  optimum  transfer  function 
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Figure 


Tracking  radar  natural  frequency,  ^ , cpt 

7.-  Effect  of  additional  filtering  in  tracking  radar;  missile- 
control  system  designed  for  optimum  Wiener  filtering. 
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E.  Missile-control  system  designed  for  optimum  Wiener 
filtering  with  additional  filtering  in  tracking  radar 

F.  Missile-conlrol  r.ystem  designed  for  both  optimum 
Wiener  filtering  and  additional  filtering 

G.  Tracking  radar  designed  for  optimum  Wiener  filtering 
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Figure  8.-  Effect  of  practical  limits  on  minimum  error;  optimum  for 

design  value. 
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Figure  9 


Tim*  constant  of  odd*d  filt*rt  second* 

- Effect  of  additional  filtering  in  tracking  radar  designed  for 
optimum  Wiener  filtering. 
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